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ABSTRACT 

Recent observations of greatly amplified magnetic fields (SB/ B ~ 100) around supernova shocks are 
consistent with the predictions of the non-linear theory of particle acceleration (NLT) , if the field is 
generated upstream of the shock by cosmic ray induced streaming instability. The high acceleration 
efficiencies and large shock modifications predicted by NLT need however to be mitigated to confront 
observations, and this is usually assumed to be accomplished by some form of turbulent heating. 
We show here that magnetic fields with the strength inferred from observations have an important 
dynamical role on the shock, and imply a shock modification substantially reduced with respect to the 
naive unmagnetized case. The effect appears as soon as the pressure in the turbulent magnetic field 
becomes comparable with the pressure of the thermal gas. The relative importance of this unavoidable 
effect and of the poorly known turbulent heating is assessed. More specifically we conclude that even 
in the cases in which turbulent heating may be of some importance, the dynamical reaction of the 
field cannot be neglected, as instead is usually done in most current calculations. 
Subject headings: acceleration of particles — shock waves — magnetic field 



1. INTRODUCTION 

The supernova remnant (SNR) paradigm for the origin 
of galactic cosmic rays is based on the assumption that 
at least ~ 10 — 20% of the kinetic energy of the expanding 
shell is converted into cosmic rays (CRs). Moreover, as 
recent observations have proved, the magnetic field (MF) 
in the s hock vicinity is amplified by a large amount (e.g. 
iBalletl (|2006| )) as would be expected if cosmic rays in- 
duce streaming instability (SI) upstream of the shock. 
We stress that such MF amplification is required to ac- 
celerate protons up to ~ 10 6 GeV. The need for a satis- 
factory and self-consistent description of these points is 
sufficient to justify the development of a NLT of particle 
acceleration. 

The develo p ments of the th e ory a r e summarized in 
(jDrurvl (|1983t l. Uones fe Ellison! (| 1991( 1 . iMalkov fc Drurvl 
(2001)). The kinetic theory for arbitrary diffusion 



coefficients () Amato fe Blasil (120051)') . and even in the 
case of self-generated MFs (lAmato fc Blasil (|2006l ) and 
IVladimirov. Ellison fc Bvkov (2006)) has been recently 
developed. 

All approaches to NL shock acceleration find that the 
large pressure of the accelerated particles decelerates 
the incoming gas, and leads to total compression fac- 
tors that scale with the Mach number of the shock as 
R t ~ MJ ~ 20 — 50. Such large shock modifications 
however are at odds with observations which are better 
fit by Rt ~ 7 — 10. The problem with larger values of Rt 
resides in both the esti mated distance betw een the for- 
ward and reverse shocks (|Warren et all (J2005D ) and in the 
fit to multifreqeuncy observations with concave spectra 
(VBK05 and references therein). The reduction in the 
compression factor is almost invariably attributed to tur- 
bulent heating (TH) in the precursor (|Berezhko fc Vblkl 
(|1997l ) and later) a s due to damping of wa ves on the 
background plasma (|McKenzie fc Vb lk (1982), hereafter 
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TABLE 1 
Parameters for 5 well known SNRs. 



SNR 


uo(km/s) 


B 2 (vG) 


P w 2 x 10 3 


Cas A 


5200 (2500) 


250-390 


32 (36) 


Kepler 


5400 (4500) 


210-340 


23 (25) 


Tycho 


4600 (3100) 


300-530 


27 (31) 


SN 1006 


2900 (3200) 


91-110 


40 (42) 


RCW86 


(800) 


75-145 


14-35 (16-42) 



Note. — For SNRs discussed by [Parizot ct al. (2003) we used 
po = 0.1rrip/cm 3 in the case of SN 1006 and po = 0.5m p /cm 3 in 
the other cases, while VBK05 provide directly P w 2- 



MKV82). In fact this mechanism was originally proposed 
in order to keep the MF amplification in the linear regime 
(e.g. SB/B <C 1), but is now commonly applied to cases 
in which SB/ B 3> 1. Unfortunately, the heating process 
is quite model dependent and even its applicability to sit- 
uations of interest for SNRs can and should be seriously 
questioned. The effectiveness of the heating process can 
easily be reduced to negligible levels or artificially ampli- 
fied to unphysical levels. 

As mentioned above, a breakthrough in the field has 
recently been provided by X-ray observations: the de- 
tection of X-ray bright filaments in the outskirts of 
some SNRs allows one to infer the strength of the MF 
in these filaments, found to be B ~ 100 — 500/iG. 
Such strong fields are generally attributed to the SI in- 
duced by CRs efficiently accelerated at the shock front, 
although alternative explan ations have been proposed 
(|Giacalone fc Jokipiil (|2007fl ). In Tab. [T] we list some 
SNRs with estimated MFs: uq is the shock velocity, B^ 
is the MF downstream of the shock as inferred from the 
X-ray brightness profile and P W 2 — B\ j '(8tt PqUq) is the 
downstream magnetic pressure normalize d to the bulk 
one. T he values of the parameters are from lParizot et all 
(2006]) and from VBK05 (in parenthesis in Tab. [J). 

We show below that for the field stregths inferred for 
SNRs, the magnetic pressure is comparable or even in 
excess of the thermal pressure of the background plasma 
and that whenever this happens the dynamical reaction 



of the field on the fluid is such that the compression 
factors are substantially reduced and fall in the range 
suggested by observations. It is crucial to keep in mind 
that, contrary to the TH, which can be either suppressed 
or amplified by changing parameters on which there is 
little or no control, the feedback of the self-generated 
turbulent MF on the plasma is not model dependent and 
must be included. 

2. DYNAMICS OF A MAGNETIZED CR MODIFIED SHOCK 

The reaction of accelerated particles upstream of the 
shock leads to the formation of a precursor, in which the 
fluid speed decreases while approaching the shock. One 
can describe this effect by introducing two compression 
factors Rt = Uq/u% and R s = U1/U2, where u is the fluid 
speed and the indexes '0', '1' and '2' refer to quantities 
at upstream infinity, upstream and downstream of the 
subshock respectively. 

The most general equations of conservation of mass, 
momentum and energy in the stationary case for a par- 
allel shock are: 

d d 

— (pu) = 0, -Q- (pu 2 + p + p c + p, w ) = 0, (1) 



8 (\ 



Tx \2 PU 



jpu 



F„ = 



dPc ro s 

- U ^- (2) 



As usual, p, u, p and 7 stand for density, velocity, pres- 
sure and the ratio of specific heats of the gas, while p w 
and F w represent the pressure and energy flux in the form 
of Alfven waves. p c is the CR pressure. The continuity 
of the distribution function of the accelerated particles 
through the subshock implies that the CR pressure is also 
continuous {p c \ — p c i)i and that the term dp c /dx gives 
null contribution when Eq. [2] is integrated from x = CP 
to x = + . All previous equations at the subshock read 
as the usual Rankine-Hugoniot relations at a magnetized 
gaseous shock. 

In order to treat the pres ence of Alfven waves co r rectly , 
we use the approach of iVainio &: Schlickeiserl (|1999f ) 
(hereafter VS99), considering two upstream wave trains 
with helicities H c = ±1, and their respective downstream 
counterparts. If SBi is a mode of the MF perturbation, 
we write the velocity perturbation as Su = —H c 



SB 

/iirp' 



Neglecting the electric field contribution, which is of or- 
der (u/c) 2 , the magnetic pressure and the energy flux, 
which for Alfven waves is the sum of the normal compo- 
nent of Poynting vector u x 5B x SB/in and the trans- 
verse kinetic energy flux p5u 2 /2, are 



- = IE«ft 



F„ = 



Z^ 



SB 2 . 
1 u. 



—- L u i +p w u, (3) 



4tt 



H c .iVA- The upstream mag- 



8tt 

having posed w, 
netic turbulen c e typi cally shows two opposite helicities 
(|Bell fc Lucekl (|2001f 0. each of which splits into a re- 
flected and a transmitted wave crossing the subshock. 
According to VS99, the transmission and reflection co- 
efficients, in the limit M 2 A ~^> R s (large Alfvenic Mach 
number), do not depend on H c and read 



R^(R S 



(4) 



For a typical supernova shock, the Alfvenic Mach num- 
ber is Ma,i = ui/va > 100, hence in the following we 



adopt these coefficients and neglect va with respect to 
the fluid velocity in Eq. [3] For each H c we therefore 
have 5B 2 /SBi = T + R = R s and thus p w2 = p w \R 2 . 

As pointed out above, the subshock can be viewed as a 
simple shock in a magnetized gas, therefore the pressure 
jump is (VS99) 



p 2 = (7 + l)R a - (7 - 1) + (7 - 1)(R S - 1)A 
Pi 7 + l-(7-l)J2 8 

with A defined as: 

A = R s + 1 [ Pw ] _ 2R S [F w ] 
R s - 1 pi R s ~ 1 pi-ui 

Using the expressions for T and R (Eq. [4]) we get 



{R s -l) 2P -^+R 



B_ ■ B, 



(5) 



(G) 



(7) 



Pi 2npi 

Following VS99, we assume that the two opposite- 
propagating waves carry MFs B± displaced in such a 

way that _B_ • B + = 0. This is not the most general con- 
figuration, but it is nevertheless rather common since it 
occurs when: 1) there is only one wave train, 2) when 
the two fields are orthogonal, and 3) on average, when 
the relative phase between the wave trains is arbitrary. 

At this point we normalize all quantities to the ones 
at upstream infinity: U(x) = u(x)/uq, P w {x) ~ Pw/pou 2 , 

and P(x) = p(x)/pqUq = —^jp—. In the latter, we used 

the assumption of adiabatic heating in the precursor and 
the conservation of mass. 

Substituting Eq. [5l Eq. [7] and the above expression for 
P(x) in the equation for momentum conservation, the 
compression factors R s and Rt are related through the 
equation 



H? 



7+1 



M 2 R] 



7+l-i? s ( 7 -l) 



1+Az 



(8) 



which is the same as the standard relation apart from 
the factor (1 + Ab), with 



A B = W [1 + R s (2/ 7 - 1)] 



(9) 



and W — P w \jP\. It is clear that the net effect of the 
magnetic turbulence is to make the fluid less compress- 
ible, noticeably reducing R t if W — P w i/P\ is of order 
1. Moreover, the pressure and temperature jumps at the 
subshock are enhanced (Eq. [5|) . 

We should notice that if one naively assumed that 
F w = 3up w everywhere, neglecting the T and R coef- 
ficients needed to satisfy Maxwell equations at the sub- 
shock, one would obtain A' = [(R s - l) 2 - 2R S ]W < 
A , leading to an incorrect pressure jump. This a 
proach, adopted bv lVladimirov. Ellison &: Bvkovi <|200 
also leads to a less marked decrease of R t , since A^ = 
VF[l + i? s (3/ 7 -2)]. 

3. CONFRONTING OBSERVATIONS 

Here we show that the magnetization levels estimated 
in SNRs as reported in Tab. [T]imply that W > 1, so that 
the dynamical feedback of the amplified MF needs to be 
taken into account. 

In Fig. Q] we plot R t versus R s for M = 100: the 
three shadowed regions represent the relation between 



J 25- 




FlG. 1. — R s -Rt for W = 1,3, 10 (intersections of solid lines with 
corresponding shadowed regions) and for \ = 200, 500 (dash-dotted 
and dotted lines). The dashed line represents the case W = (see 
text for details). 

R s and R t for fixed P w2 G [0.02,0.04] and W = 1, 3, 10: 
P w2 /R 2 s = P w i = WP gl = W (R t /R s ) J /( 7 M 2 ) . 

The three solid lines represent the relation R t — R s for 
the three given values of W as given by Eq. the dashed 
line refers to W = 0, when p w is not included. 

The compression factor lies at the intersection between 
the curve and the shadowed region for a given value of W . 
If W < 0.7 there are no intersections. This implies that 
the values of the magnetic pressure inferred from obser- 
vations require substantial MF amplification upstream, 
and that the conservation equations are affected by the 
dynamical reaction of the field. Only values W > 3 are 
compatible with the whole range 0.02 < P W 2 < 0.04 in- 
ferred from observations. This means that in order to 
account for the inferred values of B 2 , the magnetic pres- 
sure must be at least comparable to the gas pressure, and 
thus its dynamical role cannot be neglected. From Fig. 
Q] one also sees that the magnetic reaction leads to val- 
ues of Rt lower by roughly a factor ~ 2 compared with 
the case W = 0. We will comment further on this point 
below. 

Up until this point we never used the physically crucial 
point that the observed fields may be generated through 
a cosmic ray induced SI upstream of th e shock. Th e 
instability may operate in the reso nant (|Belll (|1978ah ) 
and in the non-resonant ([Belli (|2004D ) regime. 

The growth rates of these different modes can be eas- 
ily estimated only in the context of quasi-linear theory. 
Given the difficulty in deriving this information in the 
general non linear case, here we assume the following 
general relation between the pressures of CRs and wave 
pressure upstream of the subshock: 

ii = xPwi, (io) 

where £(x) = p cr (x)/poUQ is the normalized CR pressure. 
For resonant SI, one has x — Ma — uq/va, while for non- 
resonant modes, x ~ 4c/ wo- In both cases, for typical 
values of the parameters, one obtains 200 < x < 500. 

From Eq. [1] applied to the precursor, namely between 
upstream of the subshock (x = 0~) and upstream infin- 
ity, we can write: 



Rt 



1 



yM£ 



^ -1 



The physical values of R s and Rt for a given x are ob- 
tained by determining the intersection of the correspond- 
ing curve with that obtained for a given value of W at 
the subshock. Whether the solution reproduces the esti- 
mated value of P W 2 depends on whether the intersection 
falls within or outside the shadowed region for the same 
W in Fig. Q] The dash-dotted and the dotted line show 
the results for \ = 200 and 500 respectively: it is evi- 
dent that the chosen values of x allow for a consistent 
explanation of the downstream magnetic pressures as in- 
ferred from observations, and, equally important, lead 
to compression factors which ar e much lower than those 
predicted by the standard NLT (jBerezhko fc Volkl (|1997l ) 
and papers that followed). 

4. HEATING IN THE PRECURSOR 

The strong shock modification predicted by NLT when 
the magnetic pressure is ignored is usually assumed to 
be somewhat mitigated by heating o f the precursor as 
a resu lt of damping of Alfven waves (jVolk fc McKenzia 
dl98lD, hereafter VMK81, and MKV82) on the back- 
ground gas. Other phenomena (for instance acoustic 
instability) may also lead to heating of the precursor. 
In the original description, that remained basically un- 
changed to the present time, VMK81 assumed that the 
rate of damping (T) equals the rate of growth (<r) of 
Alfven waves. The main implication of this assumption is 
that the growth of the waves never reaches the non-linear 
regime, which is in fact the very reason why the mech- 
anism was invoked in the 80 's. The recent observations 
prove that waves can grow to 5B/B ^> 1. It is therefore 
at least not self-consistent to apply the standard treat- 
ment for TH to situations in which MF amplification to 
the non-linear regime takes place. In a minimal attempt 
to include faster growth one may as sume that T = aa, 
with a < 1. Following MKV82 and iBerezhko fc Ellisonl 
(1999) one can then obtain a generalized relation between 
R t and R s in the form 



Rl +1 = 
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where 
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(12) 



(13) 



P wl (l + X) = l- (11) 



which becomes equivalen t to the standard Eq. 50 of 
IBerezhko fc Ellisonl (|1999l ) only for a = 1 . Now it is easy 
to check that for typical values of R s and R t Ath > As 
if a > 3VF4t^- For instance for M A ~ 10 3 and M - 100 
one requires a to be of order unity. In this case however 
it is not easy to amplify the MF to SB ^> B^. If a is 
appreciably smaller than unity, the main process for the 
smoothening of the precursor is the dynamical reaction 
of the self-generated MF. In both cases the role of TH 
can be seriously questioned. 

A deeper look at the physical processes that may re- 
sult in the heating of the precursor make the role of TH 
even more uncertain: in the original papers of VMK81 
and MKV82 the Alfven heating was considered as a re- 
sult of non-linear Landau damping in a gas in the hot 
coronal phase. The authors reached the conclusion that 
the damping is important if Uq <C 4000fcm/s(To/5 x 
10 5 K) 1 ^ 2 , where uq is the shock velocity and Tq is the 



temperature of the unshocked gas. It is all but clear 
whether for the velocities and temperatures that apply 
to the SNRs in Tab. [TJ non-linear Landau damping is 
such to lead to a ~ 1. We stress that at the same time, 
a cannot be too close to unity, otherwise TH inhibits the 
growth of SB to the observed levels. 

Other types of turbulent heating may be at work but a 
quantitative analysis of these phenomena is lacking at the 
present time. The expression for A^h is however rather 
general, in that we did not specify the mathematical form 
of the growth and damping rates. Therefore we expect 
to draw similar conclusions in terms of the parameter a. 

This section strongly suggests that, contrary to 
the common wisdom, the most likely reason for the 
smoothening of the precursor is the dynamical reaction 
of the generated MFs rather than some form of non adi- 
abatic heating in the precursor. 

5. CONCLUSIONS 

It is well known that the effect of a MF is in gen- 
eral that of reducing the plasma compressibility. We 
showed here that when applied to a cosmic ray modi- 
fied shock, 1) this finding implies that CR induced SI 
is adequate to explain the magnetization inferred from 
X-ray observations; 2) the downstream MFs imply that 
W ~ 1 — 10, so that the field becomes dynamically im- 
portant, since this happens whenever the magnetic pres- 
sure upstream becomes comparable with the gas pres- 
sure, namely when W > 0.7 — 1; 3) the dynamical reac- 
tion of the MF reduces the compression in the precursor, 
leading to smaller (larger) values of R t (R s ) in agreement 
with the values required to explain the distance between 
forward and reverse shock and the multifrcquency obser- 
vations of several SNRs; 4) this effect comes from first 
principles, though in our calculations we restricted our 
attention to the case of Alfven waves, and is not affected 
by the huge uncertainties typical of TH; 5) an efficient 
TH may smoothen the precursor if a is close to unity, 
but in this case it is likely to inhibit the growth of the 
field to SB > B . 



Although the underlying physics is well known, the 
dynamical effect of the magnetic pressure has not been 
included in the calculations o f mult ifrequency emission 
from SNRs (|Berezhko fc Vo lk (1997) and successive pa- 
pers), so that the strong modifications predicted by NLT 
had to be compensated by assuming TH. The only ex- 
ception that we are aware o f is t he recent work by 
IVladimirov. Ellison k. Bvkovl (J2006), in which the au- 
thors perform Monte Carlo simulations of the particle ac- 
celeration process including the pressure of self-generated 
MFs. In such simulations, which represent the state of 
the art in the field, however, thermal and accelerated 
particles are treated in the same way, therefore the con- 
dition W ~ 1 could not be tested. We suspect that 
for this reason the smoothening of the precursor was at- 
tributed mostly to the backreaction of the accelerated 
particles on the field through injection. This effect is 
certainly present but as we showed here by using only 
a hydrodynamical approach, the smoothening is in fact 
mainly due to the reaction of the magnetic pressure on 
the background plasma. 

The smoothening of the precursor also results in two 
important effects: 1) spectra of accelerated particles 
closer to power laws, though the concavity which is 
peculiar of NL DSA remains evident. 2) The max- 
imum momentum of accelerated particles for given 
Mach number is predict e d to b e somewhat larger (see 
iBlasi. Amato fc Capriolil (|2007t ) for a detailed discus- 
sion). Both these effects will be discussed in detail in 
a forthcoming paper. 
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